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ABSTRACT 

The nebular spectra of the broad-lined, SN 1998bw-like Type Ic SN 2002ap 
are studied by means of synthetic spectra. Two different modelling techniques 
are employed. In one technique, the SN ejecta are treated as a single zone, while 
in the other a density and abundance distribution in velocity is used from an 
explosion model. In both cases, heating caused by 7-ray and positron deposition 
is computed (in the latter case using a Monte Carlo technique to describe the 
propagation of 7-rays and positrons), as is cooling via forbidden-line emission. 
The results are compared, and although general agreement is found, the stratified 
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models are shown to reproduce the observed line profiles much more accurately 
than the single-zone model. The explosion produced ~ 0.1 M of 56 Ni. The 
distribution in velocity of the various elements is in agreement with that obtained 
from the early-time models, which indicated an ejected mass of ~ 2.5 M & with 
a kinetic energy of 4 x 10 51 erg. Nebular spectroscopy confirms that most of the 
ejected mass (~ 1.2 M Q ) was oxygen. The presence of an oxygen- rich inner core, 
combined with that of 56 Ni at high velocities as deduced from early-time models, 
suggests that the explosion was asymmetric, especially in the inner part. 

Subject headings: supernovae: general — supernovae: individual (SN2002ap) - 
nucleosynthesis — gamma rays: bursts 



Introduction 



One of the most exciting recent discoveries in astrophysics is that the nearest ob- 



served long-duration gam ma-ray bursts 



GRB980425/SN 1998bw iGalama et al 



GRB s) are associated with sup ernovae (SNe) (e.g. 



1998; GRB030329/SN 2003dh. IStanek et all 12003 



GRB031203/SN20031w. lMalesani et al.ll2004f ). Adding to the excitement, it turned out that 
these SNe are n ot at all ordinary Spectroscopically, they are of Type Ic (no H, no He, 



weak Si lines; see lFilippenkolll997l for a review of SN classification), but their distinguishing 
feature is that unlike normal SNe Ic, they show extremely broad absorption lines dominated 
by Fe, Ca, and O. Very high expansion velocities, reaching 0.1c, are easily observed at early 
times. 

Broad lines suggest a very energetic explosion, or at least a large amount of energy 
per unit mass. Estimates of the explosion kinetic energy assuming sp herical symmetry for 
SNe 1998bw, 2003dh, and 20031w are E v , v ^ (3 - 6) x 10^ erg (e.g., 



Iwamoto et al. 1998 



Nakamura et all l200ll : Mazzali et all l2003t iDeng et all 12004 : Mazzali et all 120061 ) . This is 
more than one order of magnitude larger than in typical core-collapse SNe, which have 
-Skin ~ 10 51 erg. Ac cordingly, these exc eptionally powerful SNe have sometimes been called 
"hypernovae" (e.g., Nomoto et al .1120051 ) . In the presence of significant deviations from spher- 
ical symmetry, the se estimates may be reduced, but even then the energies remain large (e.g., 



Maeda et al.l 120061 ). The mechanism behind these events is th ought to be the collapse of th e 
stripped core of a massive star (~ 25 — 60 M ) to a black hole (IMacFadyen fe Woosleylll999l ). 



Addit ionally, evidence has been obtained that X-ray fl ashes (XRFs), the soft analogues 



of GRBs (IHeise et al. 



2006: Pian et al 



20011 ). are also linked to SNe (e.g., iFynbo et al.ll2003l : iModjaz et al. 



20061 ) . The SNe associated with XRFs also seem to be overenergetic, 
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but not as much so as GRB-SNe (jTominaga et al.l 120031 ; iMazzali et al.ll2006l ). They have 
been suggested to be the result of the collapse of the stripped carbon-oxygen core of stars of 
originally ~ 20 M & to a neutron star (IMazzali et al.ll2006l ; iMaeda et al.ll2007l). If the neutron 
star is born highly magnetized and rapidly spinning — a "magnet ar" (IDuncan &: Thompson 



19921 ; iThompson. Chang, fc Quataert 



and produce the XRF (IMazzali et al. 



2006 



2004j) — it may cause the explosion to be overenergetic 



Maeda et al.l 120071 ). 



One difference between SNe in GRBs and XRFs may be the degree of asphericity. 
GRB-SNe are thought to be highly aspherical, as GRB s are generally believed to be highly 



asymmetric phenomena (see IWoosley fc Blooml 120061 for a recent re yiew of the SN-GR B 



connection). For XRF-SNe, the degree of asphericity may be smaller (IMazzali et al.l 120071 ) 



Evidence for asphericity is not easy to glean from early-time data. The significant mixing 
outward of 5 6 Ni required to rep roduce the early rise of the light curve is a general feature of 
hypernovae (IMaeda et al.ll2003l ) and suggests an aspherical explosion, as does the connection 
of these SNe with GRBs. The best time to look for signatures of asphericity, however, is 
starting a few months after the explosion, when the SN becomes nebular, exposing the 
deepest parts of the ejecta. In that phase, both the spectrum and the light curve have a 
characteristic behavior. The light curve of SN 1998bw showed a phase o f strictly exponential 
decline, but at a rate steeper than that of 56 Co (IPatat et al.l 1200 ll ). In this phase the 
luminosity is larger than what is predicted by spherically s ymme tric models that fit the peak 
of the light curve (INakamura et al.ll200ll ). IMaeda et al.l (120031 ) showed that this behavior, 
which is actually not unusual in SNe Ic, can be explained using a modified spherically 
symmetric model where a relatively massive but slowly expanding inner core is placed in 
the center of the expanding ejecta. Although this is technically still a spherically symmetric 
model, one-dimensional explosion simulations do not predict such a density distribution, 
which su ggests that we are observ i ng an aspherical explosion, similar perhaps to the collapsar 
model of iMacFadven fc Wooslevl fjl999h . 



The nebular spectra provide more direct evidence. In SN 19 98bw the [Fe 11] lines are 
broader than the [O 1] AA6300, 6363 doublet (IMazzali et al.ll200ll ). This is also something 
that cannot be explained in the context of a spherically symmetric model of the collapse and 
explosion of a CO core, as in such a model 56 Ni (which decays into 56 Co and then into 56 Fe) 
is synthesized near the compact remnant and therefore is always located at sma ller velocities 



than O, which is left from the progenitor in the unprocessed outer layers. IMaeda et al. 



(120021 ) showed that this configuration can be obtained in an axisymmetric explosion. In 
such an explosion, most of the E^ iQ is released along the "jet" axis, which is probably linked 
to the launching of the GRB, and 56 Ni is synthesized along that direction and ejected at a 
high velocity. Away from the jet axis, however, less kinetic energy is deposited and burning 
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proceeds much less efficiently. Therefore, in these directions large amounts of unburned 
are ejected at low velocity. For a near-polar viewing angle, this scenario naturally leads to 
broad [Fe n] lines and a narrow [O i] line at late times. For an equatori al view, on the othe r 
hand, the [O i] line should have a characterist ic double-peaked pr ofile (IMaeda et al.l 120021 ). 
as was indeed observed in the SN Ic 2003jd (IMazzali et al.ll2005l ). Fur ther evidence that 



SNe Ic are aspherical also comes from their high degree of polarization (IWang et al.l 12001 



Leonard et al.ll2002l ; iFilippenko fe Leonard! 120041 ; iLeonard fe Filippenkoll2005l ). 



For XRF/SNe, the case for aspherici ty is much weak er. The early-time spectra of 

SN2006aj did not show very broad features (IPian et al. hood ). Radio observations (ISoderberg et al 
20061 ) suggest a very broad opening angl e (> 60°) . The broa d nebular lines are not inconsis- 
tent with spherically symmetric ejecta (IMazzali et al.l 120071 ) . Nevertheless, the presence of 
emission at velocities below 2000 km s -1 is at odds with the predictio n of one- dimensional 
explosion models, that a density hole is present at the lowest velocities (IMazzali et al.l 120071 ; 
Maeda et alll2007h . 



Besides those clearly associated with GRBs or XRFs, other energetic SNe Ic have been 
observed. Energetic SNe Ic, with or without a GRB, are recognized from the extreme width 
(up to ~ 0.1c) of their spectral lines at early times, and have sometimes also been called hy- 
pernovae (although their total energy is not always much higher than that of normal SNe Ic). 
The nearest such object ever observed was SN2002ap in M74. SN2002ap was immediately 



recog nized as a broad-lined event, similar spectroscopically to SN1998bw (IKinugasa et al. 



20021 ). and was therefore intensively observed. SN 2002ap remained much less luminous than 
SN1998bw (IMazzali et alil2002l : iGal-Yam. Ofek. fe Shemmerlbooi ). and it was not signifi- 
cantly more l uminous than the average SN Ic. Also, SN 2002ap wa s not seen in association 
with a GRB Jrlurlev et al.l I2OO2I : iGal- Yam. Ofek. fe Shemmerll2002h . 



Mazzali et al.l (120021 ) modelled the early-time light curve and spectral evolution of 
SN2002ap and derived values of the ejected mass (M e j = 2.5 — 5M ), the kinetic energy 
of the explosion (E^ m — 4 — 10 x 10 51 erg), and the mass of 56 Ni synthesised in the explo- 
sion (M( 56 Ni) = .07M^). While the mass of 56 Ni is similar to that of "normal" SNe Ic 
JSauer et al.l l2006h . and much smaller than that of hypernovae linked to GRBs, both M e - } 
and .Ekin are intermediate between those of normal SNe Ib/c and GRB-SNe. The uncertainty 
in the values of M e j and E^i n follows from assuming either the absence of any He envelope 
(lower bound) or the presence of a maximal, 2.5M helium envelope, the presence of which 
would not affect the light curve significantly but is at the same time not supported by the 
spectral appearance of SN2002ap. Inte restingly, the light-c urve behavior at advanced phases 
was very similar to that of SN 1998bw (ITomita et al.ll2005l ). suggesting a common nature for 
the two events. 
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Mazzali et al.l (120021 ) argued that SN 2002ap was the ~ 5 M & carbon-oxygen core of a 
star of initially ~ 25 M that collapsed to a black hole. On the other hand, no GRB was 
detec ted, only weak radio emission was detect ed from the SN (IBerger. Kulkarni. &: Chevalier 
2002 ). and the X-ray signal was very weak ( Soria. Pian. &: Mazzali 20031 ). suggesting that 
little or no relativistic ejecta were produced. This may be related to the relatively small 
mass of the collapsing star, which may have been close to the minimum required to form a 
black hole. It is therefore interesting to explore further the nature of this object. As in the 
case of SNe 1998bw and SN2006aj, this can be done by modelling the nebular spectra. 



Foley et al.l (120031 ) published a series of late-time spectra of SN 2002ap, while iKawabata et al 
( 120021 ) published a single spectrum. The nebular phase in SN2002ap start ed to develop 



rathe r early, as in SN 1998bw but unlike other hypernovae such as SN 1997ef (IMazzali et al. 
2004J). At an age of ~ 4 months, the SN spectra were already fully nebular (IFoley et al. 



20031 ). The optical spectra were dominated by very strong [O i] AA6300, 6363 emission 
and showed strong Ca II emission lines, similar to SN 1998bw. Unlike the case of SN 1998bw 
however, the [Fe ll]-dominated blend near 5200 A is rather weak. This is not su rprising, since 
SN 19 98bw is thought to have produced ~ 4 times as much 56 Ni as SN 2002ap (IMazzali et al. 



20061 ). Still, the [Fe n] lines in SN2002ap are sufficiently pronounced that modelling can 
be meaningfully attempted. One peculiarity of SN 2002ap is the great strength of the Mg i] 
AA4571 line, which may suggest that the envelope of SN 2002ap was more thoroughly stripped 
th an that of SN 1998bw . This would favor the lower bound of the mass and energy estimates 
of IMazzali et al.l (120021 ). which were obtained assuming the absence of any significant He 
layer. 



2. The Spectral Dataset 



In this paper, we model all the spectra published by IFoley et al.l ( 120031 ). Additionally, 
we present and model three additional late-time spectra. Two of these spectra were obtained 
with the Subaru telescope. The other spectrum was obtained at the National Astronomical 
Observatory of China (NAOC, formerly Beijing Astronomical Observatory, BAO). 

The Subaru spectra were obtained on 2002 Jun 7.6 (UT dates are used throughout this 
paper) and Sep 15.6 with FOCAS. For both spectra, the exposure time was 240 s, and grism 
B300 and Y47 order-cut filters were used. A 0?8 wide slit was used for the June spectrum 
and the resulting spectral resolution (A/AA) is ~ 650. For the September spectrum, a 2 / .'0 
slit was used, but the spectral resolution is still similar to that of the June spectrum because 
the seeing size was / .'7-0 / . / 8. The flux was calibrated using observations of either Feige 110 
or G191B2B obtained on the same night as the target spectra. 
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The BAO spectrum was obtained on 2002 July 11.8 with the NAOC 2.16 m telescope 
at Xinlong Observatory (then BAO). The observations were carried out with an OMR (Op- 
tomechanics Research, Inc.) spectrograph, using a Tektronix 1024 x 1024 pixel CCD as 
the detector. A grating of 300 g mm -1 was used, which provided a spectral resolution of 
~ 10 — 11 A. The exposure time was ~ 40 min. 

The newly published spectra are shown in Figure 1. The NAOC spe ctrum is rea s onabl y 
consistent with the almost contemporaneous Lick spectrum published by iFoley et al.l (120031 ) , 
and therefore it was modelled but is not explicitly shown in the Figures in the following 
sections. 

The 11 spectra we used for modelling cover a time span of almost 9 months, from June 
2002 through February 2003, corresponding to SN ages of 4-13 months given that the SN 
exploded on 29 January 2002 (IMazzali et al.ll2002l ). Unfortunately, the wavelength coverage 
is not uniform. In particular, the Subaru spectra extend over a shorter range (4700-9000 A) 
than the other spectra, missing the important Mg i] line. Models based on the Subaru spectra 
are therefore somewhat less reliable than those based on the Lic k spectra. All spectra were 
calibrated by using the available photometry (IFoley et al.l 120031 ) with the exception of the 
last three, for which an extrapolation of the light curve was used. 



3. Modelling Technique 



In order to model the nebular s pectra of SN2002ap , we used our non-local thermody- 
namic equilibrium (non-LTE) code (IMazzali et al.l 120011 ) . The code computes the heating 
of the gas following the deposition of 7-rays and positrons emitted by the decay of 56 Co 
into 56 Fe. Heating is balanced by cooling via nebular line emission. The emission rate in 
each line is comp uted solving a non-LTE matr ix of level populations (lAxelrodlll980l ). In the 
original version (IRuiz-Lapuente fc Lucy! |l992j) , a homologously expanding nebula of finite 
extent, uniform density, and uniform composition is assumed, and the emission spectrum is 
obtained assigning to all lines a parabolic profile, bounded by the velocity of the outer edge 
of the nebula. 



Together with this classical version, a more advanced — although more model-dependent 
- version has been developed and is used in this work. In this new version, stratification 
in density and abundance is adopted. The density profile is taken from explosion models, 
and 7-rays and positrons are emitted at various depths according to the distribution of 56 Ni. 
Their propa gation and deposi t ion is followed using a Monte Carlo scheme similar to that 
discussed by ICappellaro et al.l (119971 ) for their light-curve models. A constant 7-ray opacity 



- 7- 



(/c 7 = 0.027 cm 2 g -1 ) and a constant positron opacity (« e + = 7cm 2 g~ 1 ) are assumed. The 
heating and cooling of the gas are then computed in non-LTE in each radial shell, and so is 
the line emissivity. The line profiles from each shell are assumed to be truncated parabolas, 
the inner truncation point corresponding to the inner boundary of the shell considered. These 
truncated parabolas are then summed to produce the emerging spectrum. Line profiles thus 
depend on the density and abundance distributions. Therefore, this approach constitutes a 
test of the explosion models which are used to simulate the light curves and the early-time 
spectra, especially for what concerns the innermost part, which is only visible directly in the 
nebular phase. 

Because the stratified code depends more directly on the adopted explosion model than 
the one-zone version, it is useful to compare the results of the two approaches. The next two 
sections deal with both cases in turn. 

One important ingredient for the modelling is the assumed distance and reddening 



Vinko et al. 


f 


2004 





supergiant stars, obtained a distance modulus of \i — 2 9.32 mag for M7 4 , and \x = 29.50 



Vinko et all faooAi comm ent, this value was derived using 



mag for the M74 group, confirm ing previous result s by ISohn &: Davidgd (119961 ). who pre- 
ferred to use red supergiants. As 
the Galactic absorption maps of 
absorption in the B band of Ar = 0.13 m ag 



Burstein fc Heiled (119821 ). which give a value of the total 



On the other hand, the more recent maps 

0.301 



published by ISchlegel. Finkbeiner. fc David (119981 ) predict a larger absorption, A B 
mag, leading to a reduced distance, \i = 29.15 mag. 



As for the reddening, combining the Galactic extinction of lSchlege 



2002 



(1998), E(B— V)g = 0.075 mag, with the small reddening within M74 (ITakada-Hidai. Aoki. fc Zhao 



, Finkbe iner. fc Davis 



E(B - V) H = 0.020 mag), a total value of E(B - V) = 0.09 mag is used. 



Clearly, the uncertainty associated with the distance is the largest contribution to the 
overall uncertainty on the luminosity of SN 2002ap. Therefore, in this paper we have chosen 



to adopt the same valu es used in previous modelling and analysis papers (IMazzali et al. 



2002 



Yoshii et al.l bood . (ji = 29.50 mag, E(B - V) = 0.09 mag)) in order to make the 



results immediately comparable. The effect of adopting a different value for the distance is 
discussed below. 
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One-Zone Models 



The driving parameters for the one- zone fits are the line width, the mass of 56 Ni, and the 
masses of the other elements. Fitting the widths of the complex [Fe n] blend near 5200 A gives 
a measure of the distribution of 56 Ni, since most Fe is the product of the decay of 56 Ni. This 
blend can be fitted for a nebular velocity of between 6200 and 5200 km s -1 . The line width 
slowly decreases over the period con sidered, which is in good agreement with what was found 
for SN1998bw (IMazzali et al.ll200ll ). although in that case the velocities were significantly 
larger. This behavior indicates that the outer parts of the ejecta are becoming progressively 
more difficult to excite as the density decreases owing to the expansion. Unlike the case of 
SN 1998bw, however, the velocity that is required to fit the [Fe n] lines is also appropriat e 
for the [O i] line. This is typical of lower-energy SNe Ic (e.g., SN 19941. ISauer et al.ll2006l ). 
and indicates that the degree of asphericity of the explosion is small. The only apparent 
deviation from sphericity is the presence of narrow emission spikes in the co res of [O i] 
AA6300, 6363 and Mg i] A4571, at velocities bel ow 2000 km s" ^Folev et al.l l2003l. Figure 17). 
Such a feature was also observed in SN2006aj (IMazzali et al.l 120071 ). 



The masses of the other elements are determined by fitting the various emission lines. 
One caveat is that two intermediate-mass elements that are expected to be abundant in the 
SN ejecta are only diagnosed by a single line. The only line of sulfur is [S i] AA4069, which 
is blended with various Fe lines. Silicon, a highly abundant element, also has no strong or 
isolated line at optical wavelenghts. A shoulder on the red side of [O i] AA6300, 6363 may be 
due to [Si i] A6527. We have determined the masses of Si and S by trying to reproduce those 
two features, but the uncertainty involved is large, affecting also the estimate of the mass in 
the ejecta. In fact, both silicon and sulfur have strong lines in the infrared ([Si i] 1.61, 1.65 
/im, [S i] 1.08, 1.13 /nil), which act as efficient coolants. Thus, increasing the mass of these 
two elements requires all other masses, including that of 56 Ni, to be increased. Additionally, 
since in the later spectra it is not possible to reproduce accurately the red shoulder of the 
[O i] line, there may be some doubt as to the actual contribution of silicon to that feature. 
More accurate statements about the masses of these elements would require the availability 
of infrared spectra. 

As in SN1998bw, the spectra of SN2002ap show lines of [Fe n], but not [Fe ill]. This 
implies a rather low degree of ionizati on, which can be rep roduced assuming significant 
clumping to favor recombination. As in IMazzali et al.l (120011 ). a volume filling factor of 0.1 
was adopted in all models to achieve this. The models with stratified density and composition 
discussed in the next section are useful to address this issue. 



The results are summarized in Table 2, and the time series of the spectral fits is shown 
in Figures 1 and 2. The electron temperature and density in the nebula decrease with 
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time. The inferred value of M( 56 Ni) is ~ 0.10 M & . Small oscillations around this value 
in the various fits are probably due to a combination of inconsistent flux calibration and 
incomplete wavelength coverage, in particular for spectra that do not include the Mg i] line. 
Only in the two very late-time spectra does the 56 Ni mass increase significantly, but these 
spectra are calibrated in flux using an extrapolation of the observed photometry, making the 
results less reliable. 

In order to fit the spectra, we require an ejected mass of ~ 1.6 — 1.9 M & below a velocity 
of 5500 km s" 1 . Of t his mass, ~ 0.7 M is oxygen. The density profile (CO100/4) used by 
Mazzali et al.l (120021 ) to fit the early-time light curve and spectra only contained ~ 0.8 M 



be low the same veloci ty. Interestingly, the larger mass is in good agreement with the results 
of iMaeda et al.l (120031 ). who used a two-component, one-dimensional density distribution to 
reproduce the light curve of SN2002ap and obtained a 56 Ni ma ss of 0.08 M and an ejected 
mass of 1.6 M below 5750 km s _1 (see also lTomita et al.ll2005l ). Although the mass at low 
velocity may be somewhat overestimated by the one-zone models, which overfit the flux at 
the lowest velocities where the observed [O i] line has a narrow core, the need for additional 
mass at low velocity is clear. Simila r results were also obtained for other SNe Ic (e.g., 
Sauer et al.l 120061 ; iMazzali et al.l 120071 ) . They suggest the presence of a dense, oxygen-rich 
core in the ejecta, which would be most naturally explained as the result of an aspherical 
explosion. 

One place where our models fail is the emission near 7700 A. This is likely to be O I 
A7774. The excitation temperature of this line is significantly higher than the nebular tem- 
peratures obtained in our models, suggesting that the flux in the line may be caused by 
recombination or non-thermal excitation by fast electrons. 

Therefore, while the results f rom the one-zo n e mod els generally confirm findings from the 
light curve an d spectral study of Mazzali et al.l (120021 ). they suggest that the improvements 



introduced by lMaeda et al.l (120031 ) are realistic, and that a simple one-dimensional explosion 



model is not sufficient to explain the behavior of SN 2002ap. 



5. Multi-Zone Models 



We have modelled the same 11 spectra of SN2002ap with the multi-zone code de- 
scribed in Sec t ion 2. As a model of the explosion, we selected CO100/4, which was used by 
Mazzali et al.l (120021 ) to fit the early-time light curve and spectra of SN 2002ap. 



Our multi-zone code is still one-dimensional, but by allowing us to investigate the el- 
ement distribution in velocity space at an epoch when the ejecta are fully transparent, it 
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can give us indirect clues of possible asymmetries through unusual abundance distributions. 
Therefore, although we started our models using the density profile and the abundances of 
CO 100/4, we allowed the abundances to vary so as to achieve the best possible fit to the 
data. Repeating this exercise over several epochs ensures the validity of the results. 

Also, because the ejecta are assumed to be fully transparent, once the abundance dis- 
tribution has been defined from fitting one spectrum, it should not need to be changed to 
reproduce the other epochs. Ideally, this should be done with the earliest nebular spectrum, 
since this is likely to show emission from a broader range of velocities than later spectra, as 
discussed in the previous section. 

Since the first spectrum of our series does not cover the Mg i] A4571 line, which is an 
important coolant, we used as reference the Lick Observatory 8 June 2002 spectrum, which 
offers very broad wavelength coverage at a high signal-to-noise ratio. The model used to fit 
that spectrum gave good results for the other spectra, with only small changes in composition, 
which may be attributed to the variable observing conditions and flux calibration. 

The series of our synthetic spectra is shown in Figures 4 and 5. In order to reproduce 
the width of the [Fe n] lines, the mass fraction of 56 Ni must be small (~ 1%) at 9000 < v < 
17, 000 km s _1 , and then slowly increase inward, reaching ~ 10% at 3000 < v < 6000 km s _1 . 
At the lowest velocities, below 3000 km s _1 , the abundance drops again to ~ 2%. For the 
highest velocities (v > 17,000 km s _1 ), the value of the abundance derived from the early- 
time modelling (~ 10%) was used, but this does not affect the nebular spectra. The oxygen 
mass fraction, on the other hand, does not change significantly with velocity. It is ~ 65% 
at v > 10,000 km s™ 1 , decreasing to ~ 40% at 2000 < v < 7000 km s" 1 , a region that is 
dominated by 56 Ni, Si, and S. The abundance of S was determined from the weak [S i] A4069 
line, and the Si abundance was assumed to have a constant ratio of 3:1 with respect to 
that of S. With this assumption, the emission on the red side of the [O i] line is not fully 
reproduced, suggesting that [Si n] is not the only line contributing to it, or that a larger 
mass, or a larger Si/S ratio should be adopted. Other elements behave as expected. The Mg 
and Ca abundances also increase inward, while C decreases with O. 

One problem using model CO100/4 is that the strong, narrow cores of [O i] AA6300, 6363 
and Mg i] A4571 cannot be reproduced. The characteristic velocities of these cores are in fact 
< 2000 km s" 1 , but CO100/4 has a density "hole" at v < 3000 km s~ l . Therefore, following 



Maeda et al.l (120031 ). we added a dense inner region at v < 3000 km s _1 . The composition of 
this region must be dominated by O, with C, Mg, and other intermediate-mass elements also 
present so that the line cores can be reproduced. However, only little 56 Fe must be present 
there, or the [Fe il] feature would be negatively affected. This inner zone contains only 
~ 0.2 M Q of material, mostly oxygen. The total 56 Ni mass below 3000 km s -1 is ~ 0.01 M & , 
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in agreement with the result of iMaeda et al.l (120031 ). The oxygen mass below 5750 km s 1 in 



these models is 0.35 M & , which is significantly less than in the one-zone models. The total 
56 Ni mass is, however, slightly larger (~ 0.11 M Q ). The total ejected mass is M e j « 2.5 M , 
of which ~ 1.3 M & is oxygen. 



The main results obtained from the multi-zone models are recapped in Table 2. The 
last two spectra again require a somewhat larger 56 Ni mass. Because of the uncertain flux 
calibration of these spectra, we do not regard this as an inconsistency of our modelling 
procedure. The oxygen mass is, however, much more stable than in the one-zone models, 
showing that the multi-zone approach is more accurate, and that the explosion model that 
we used is a good representation of the ejecta of SN 2002ap in the late phase as well as the 
early phase. 

One further point concerns clumping. Introducing a density gradient should reduce 
the ionization degree at low velocity. Thus, it may be expected that in the stratified models 
clumping would not be necessary. However, even in these models, if clumping is not used the 
Fe ionization is too high at all velocities. A filling factor of 0.1 was therefore used in all zones 
with v > 2000 km s" 1 — that is, wherever 56 Ni is present. No clumping is necessary in the 
O-dominated inner zone, below 2000 km s _1 . There are a number of possible explanations 
for this. One is that clumping does indeed exist. The lack of direct evidence (e.g., narrow 
emission features at different velocities) suggests that any clumps must be small and evenly 
distributed in velocity space. Another possibility is that 56 Ni is located mostly near the 
jet direction in an axisymmetric explosion. In this case, the actual density of Fe in the 
effectively emitting volume in the nebular phase would be higher than what is suggested 
by a spherically symmetric approach like the one used here. Clumping would then be a 
proxy for this aspherical distribution. The total mass estimate would not be affected, since 
in the nebular phase each emitting ion contributes to the line profile, independent of its 
position. An axisymmetric configuration was indeed suggested for the outer part of the 
ejecta as a possible interpretati on of the detection of significant line pola rization in the 
early-time spectra of SN2002ap ( Kawabata et al. 2002 ; Leonard et al.l I2OO2I ) . If significant 



asymmetry affected the inner part of the ejecta, the profiles of th e nebular em i ssion lines 



would depend sensitively on viewing angle, as demonstrated by IMaeda et al.l (120021 ) for 
SN 1998bw. However, in this case the lack of a coincident GRB and of dramatic signatures 
of asphericity in the nebular spectra suggest that any deep asphericity was most likely weak. 
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6. The Geometry of the Ejecta of SN 2002ap 



We have modelled a series of nebular spectra of the broad-lined SN2002ap, using a 
one-zone non-LTE code and a multi-zone code where 7-ray deposition and line emission are 
computed more accurately based on a pre-defined explosion model. Both simulations yield 
a good reproduction of the [Fe 11] lines, but the narrow cores observed in [O 1] AA6300, 6363 
can only be reproduced assuming the presence of a dense, O-dominated inner zone. In the 
one-zone synthetic spectra, a mass excess of ~ 0.7 M is needed to reproduce qualitatively 
the [O 1] emission line. In the multi-zone model, which reproduces the line profiles much more 
accurately, the mass excess is smaller, ~ O.2M , but the presence of oxygen-dominated ma- 
terial moving at low velocity is confirmed. This material is not predicted by one- dimensional 
explosio n models, but its p resence was suggested in various hypernovae based on light-curve 
studies f lMaeda et aDl2003h . 

Similar results were al so obtained for the normal SNIc 19941 ( Sauer et al. 20061 ) as well 
as for the XRF/SN 2006aj (IMazzali et a l]2 QQj). SN 19941 shows em ission in the [O 1] AA6300, 
6363 line at velocities below 2000km s _1 (jFilippenko et al" 1995 ). while the corresponding 
one-dimensional explosion model predicts the absence of material at low velocities. The 
rounded profile suggests a smooth distribution of oxygen, which does not by itself indicate 
an aspherical explosion, but the very presence of material does, and it can be modelled adding 
O.2M of oxygen-dominated mass at low velocity. SN2006aj shows a peak in [O 1] AA6300, 
6363 at velocities below 2000 km s -1 , although not quite as sharp as in SN2002ap (Figure 
4). Thi s was modelled for a mass excess of ~ O.7M over the one-dimensional explosion 



model (IMazzali et al.l 120071 ). The presence of a dense inner zone in all SNe Ic suggests that 
these explosions are intrinsically aspherical, as is also deduced from polarization studies 
(I Wang et al.l l200ll ; iLeonard et al.l 120021 ; iFilippenko fe Leonard! 12004 ; iLeonard fe Filippenko 
2005 



Figure 6 shows a comparison of the [O 1] AA6300, 6363 line in SN2002ap, and in 
SNe2006aj, 1998bw, and 2003jd (Mazzali et al. 2007, 2001, and 2005, respectively). The 
line in SN2006aj has a full-width at half- maximum intensity (FWHM) of ~ 8000 km s _1 , 
while for the other three SNe the FWHM is ~ 6000 km s _1 . The line profile in SN2006aj is 
not very sharp, indicati ng; that the nebula is to a good approximation spherically symmet- 
ric (IMazzali et al.l 120071 ). In this case the [O 1] line velocity should correspond to the real 
expansion velocity. The sharp profile in SN 1998bw was interpreted as a disk-like distribu- 
tion of oxygen viewed from a near-polar direction. Although the low observed velocity is 
partly the result of a projection effect, the actual velocity of the oxygen (~ 8000 km s" 1 ) is 
smaller than the velocity of Fe JMaeda et alJbooj ). The [O 1] line of SN2003jd has a width 
similar to that of SN 1998bw, but the double-peaked profile suggests that we are viewing 
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the oxygen-rich disk from near its plane (IMazzali et al.ll2005l ). The case of SN2002ap seems 
intermediate between those of SNe 1998bw and 2006aj. The line can be separated into two 
components. A broad base, of width similar to that in SNe 1998bw and 2003jd, is present 
down to velocities of ~ 3000 km s -1 and is well fitted by a one- zone model. Its profile is 
reminiscent of that of SN2003jd, but the presence of the sharp core makes it impossible to 
distinguish between a flat-topped profile, characteristic of a shell-like spherical distribution, 
and a double-peaked profile as in a disk viewed edge-on. This core is narrower than even 
the corresponding profile of SN 1998bw. 

There may be two explanations for this type of profile. One is that in SN2002ap the 
outer part of the oxygen, from v ~ 3000 out to v ~ 8000 km s _1 , has a distribution that 
is not far from spherical, giving rise to a flat-topped profile. Below ~ 3000 km s _1 , on the 
other hand, we may be seeing a situation qualitatively similar to that of SN 1998bw, with 
the oxygen distributed in a disk-like region oriented not far from face-on. The size and mass 
of this disk-like region are, however, much smaller in SN2002ap than in SN1998bw. The 
other possibility is that the outer part reflects a disk-like distribution of most of the oxygen 
viewed not far from its plane, like in SN2003jd. The broad component of the profile would 
then be double-peaked rather than flat-topped, and the inner component would contain more 
mass, so that its emission could fill the valley between the two peaks. The sharpness of this 
inner-component emission may again indicate a disk-like distribution viewed face-on, but in 
this case two main orientations would exist in the ejecta, which may not be easy to explain 
in terms of the collapse mechanism. Alternatively, the inner region could simply be sharply 
peaked in density, as reproduced by our multi-zone model, but not grossly aspherical. 

In both scenarios the lack of a GRB may be justified. If the inner explosion was as- 
pherical, but the asphericity did not propagate to the intermediate part of the ejecta (3000 
to 8000km s _1 ), it is unlikely that a beam of material moving at relativistic velocities was 
produced, yet the overenergetic explosion could have caused enough material to move at high 
velocities (~ 30000 km s _1 ), giving rise to the broad absorption features in the early spectra. 
The dis tribution of this material may not have b een very far from spherical: polarization 
studies (jrCawabata et al. 2002 ; Leonard et al. 2002 ) suggest a degree of asphericity of ~ 10%. 
On the other hand, if the outer explosion was aspherical, a relativistic jet may have been 
produced, but it would not have been pointi ng to us. This latter possibili ty is, however, not 
supported by the la ck of an X-ray detection (ISoria^ Pian. fc Mazzalill2003l ) and the weakness 
of the radio signal (IBerger. Kulkarni. fc Chevalierl 120021 ). so we tend to favor an explosion 
that was not far from spherical, except in the innermost part. 

Both simulations agree that the 56 Ni mass ejected by SN2002ap is y 0.10 M^, a slightl y 
larger value than that predicted by studies of the early-time light curve (IMazzali et al.ll2002l ). 
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but in agreement with the late-time light curve jMaeda et alil2003l : iTomita et al-lbood ). The 
multi-zone model confirms that only a small fraction of this (~ 0.01 M & ) is located in the 
innermost zone, a clear indication that the explosion was asymmetric to some degree. 

Other global values that can be derived rather accurately from a one- dimensional study 
of the nebular spectra (M( 56 Ni) « 0.11 M , M e j « 2.5 M ) are small compared to those of 
other broad-lined SNe, confirming previous results based on the early-time light curve and 
spectra. These values point to a progenitor less massive than for objects such as SN 1998bw 
or SN 2003dh, although still a massive star (M w 22 — 25 M ) that could have produced 
a black hole. This also suggests that the "hypernova" properties of SN2002ap were less 
extreme, including asymmetry and the launching of a GRB. 

The estimated mass of 56 Ni would be smaller if we adopted a shorter distance. If the 
distance modulus was fi = 29.15 mag, the 56 Ni mass would be reduced by ~ 15%. The mass 
and energy of the ejecta would also be similarly reduced, but the change is sufficiently small 
that the shape of the light curve would not be affected. 

Indirect indications of asymmetry can be obtained in our one- dimensional study, but 
the details of the exact properties are not. The development of axisymmetric explosion 
models suitable for SN2002ap, and the study of th e nebular spectra using three-dimensional 
spectrum synthesis codes as in iMaeda et al.l (120021 ). are currently planned. 



The work of A.V.F.'s group at UC Berkeley is supported by National Science Foundation 
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Fig. I. — The nebular spectra of SN2002ap obtained at Subaru and NAOC. The two earlier 
spectra have been shifted upward by adding a constant. The zero level of the flux is shown 
by a thin horizontal line. 
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Fig. 2. — The nebular spectra of SN2002ap obtained from 2002 June to 2002 September, 
compared to the synthetic spectra computed with the one- zone model (dotted/red). Except 
for the lowest one, all spectra have been shifted upward by arbitrary amounts. Long tick 
marks on the left ordinate axis show the zero flux level for each spectrum. [See the electronic 
edition for a color version of this figure] 
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Fig. 3. — The nebular spectra of SN 2002ap obtained from 2002 October to 2003 February, 
compared to the synthetic spectra computed with the one- zone model (dotted/red). Except 
for the lowest one, all spectra have been shifted upward by arbitrary amounts. Both the 
lowest spectrum and the corresponding model have been multiplied by a factor of 2 for 
display purposes. Long tick marks on the left ordinate axis show the zero flux level for each 
spectrum. [See the electronic edition for a color version of this figure.} 
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Fig. 4. — The nebular spectra of SN2002ap obtained from 2002 June to 2002 September, 
compared to the synthetic spectra obtained with the multi-shell model (dotted/red). Except 
for the lowest one, all spectra have been shifted upward by arbitrary amounts. Long tick 
marks on the left ordinate axis show the zero flux level for each spectrum. [See the electronic 
edition for a color version of this figure] 
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Fig. 5. — The nebular spectra of SN 2002ap obtained from 2002 October to 2003 February, 
compared to the synthetic spectra computed with the multi-shell model (dotted/red). Except 
for the lowest one, all spectra have been shifted upward by arbitrary amounts. Both the 
lowest spectrum and the corresponding model have been multiplied by a factor of 2 for 
display purposes. Long tick marks on the left ordinate axis show the zero flux level for each 
spectrum. [See the electronic edition for a color version of this figure.} 
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Fig. 6. — A comparison of the [O i] AA6300, 6363 line of SN e 2002ap (so l id gre y/black line), 
2006aj (long dashed/blu e line), 1998bw (do tted/red line, IPatat et all boOlh . and 2003jd 
(short dashed/green line. iMazzali et al.ll2005l ). Spectra have been scaled by various amounts 
to facilitate comparison. [See the electronic edition for a color version of this figure.} 
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Table 1. Parameters of the One-Zone Synthetic Spectra 



UT Date 


SN epoch 


Telescope 


velocity 


M( 56 Ni) 


M(O) 


M(Si) 


M(tot) 




log (n e / 




[days] a 




[km s ] 


[M e ] 


[Mq] 


[Mq] 


[M ] 


[K] 


(g cm -3 )) 


7 June 2002 


129 


Subaru 


6200 


0.090 


0.71 


0.60 


1.78 


5511 


7.58 


8 June 2002 


130 


Lick 3 m 


6200 


0.092 


0.62 


0.50 


1.58 


5488 


7.56 


11 July 2002 


163 


Lick 3 m 


5750 


0.088 


0.71 


0.57 


1.73 


5051 


7.34 


11 July 2002 


163 


NAOC 


5750 


0.090 


0.71 


0.58 


1.73 


5082 


7.35 


9 Aug 2002 


192 


Lick 3 m 


5550 


0.093 


0.66 


0.60 


1.69 


4751 


7.17 


15 Sept 2002 


229 


Subaru 


5550 


0.098 


0.70 


0.74 


1.90 


4360 


6.95 


1 Oct 2002 


245 


Lick 3 m 


5550 


0.102 


0.67 


0.78 


1.86 


4209 


6.86 


8 Oct 2002 


253 


Lick 3 m 


5550 


0.103 


0.66 


0.70 


1.88 


4160 


6.82 


6 Nov 2002 


281 


Lick 3 m 


5550 


0.110 


0.71 


0.80 


2.04 


3903 


6.69 


7 Jan 2003 


343 


Lick 3 m 


5450 


0.145 


0.61 


0.80 


1.93 


3422 


6.49 


27 Feb 2003 


394 


Lick 3 m 


5200 


0.170 


0.62 


0.80 


1.92 


3092 


6.37 



a The epoch is given from the putative date of explosion, 29 Jan 2002 llMazzali et al.ll2002h . 



Table 2. Parameters of the Multi-Zone Synthetic Spectra 



UT Date 


SN epoch 


Telescope 


Af( 56 Ni) 


M(O) 


M(Si) 




[days] a 




[M Q ] 


[Mq] 


[M ] 


7 June 2002 


129 


Subaru 


0.11 


1.21 


0.44 


8 June 2002 


130 


Lick 3 m 


0.11 


1.13 


0.58 


11 July 2002 


163 


Lick 3 m 


0.11 


1.33 


0.31 


11 July 2002 


163 


NAOC 


0.11 


1.33 


0.31 


9 Aug 2002 


192 


Lick 3 m 


0.12 


1.28 


0.36 


15 Sept 2002 


229 


Subaru 


0.11 


1.34 


0.29 


1 Oct 2002 


245 


Lick 3 m 


0.11 


1.29 


0.36 


8 Oct 2002 


253 


Lick 3 m 


0.11 


1.29 


0.33 


6 Nov 2002 


281 


Lick 3 m 


0.11 


1.29 


0.31 


7 Jan 2003 


343 


Lick 3 m 


0.13 


1.28 


0.36 


27 Feb 2003 


394 


Lick 3 m 


0.13 


1.35 


0.30 



a The epoch is giv en from the putative date of explosion, 29 Jan 2002 
l lMazzali et al.ll2002ll . 



